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The idea that dark matter forms part of a larger dark sector is very intriguing, given 
the high degree of complexity of the visible sector. In this paper, we discuss lepton jets as 
a promising signature of an extended dark sector. As a simple toy model, we consider an 
O(GeV) DM fermion coupled to a new U (1)' gauge boson (dark photon) with a mass of order 
GeV and kinetically mixed with the Standard Model photon. Dark matter production at the 
LHC in this model is typically accompanied by collinear radiation of dark photons whose 
decay products can form lepton jets. We analyze the dynamics of collinear dark photon 
emission both analytically and numerically. In particular, we derive the dark photon energy 
spectrum using recursive analytic expressions, using Monte Carlo simulations in Pythia, and 
using an inverse Mellin transform to obtain the spectrum from its moments. In the second 
part of the paper, we simulate the expected lepton jet signatures from radiating dark matter 
at the LHC, carefully taking into account the various dark photon decay modes and allowing 
for both prompt and displaced decays. Using these simulations, we recast two existing 
ATLAS lepton jet searches to significantly restrict the parameter space of extended dark 
sector models, and we compute the expected sensitivity of future LHC searches. 


I. INTRODUCTION 

For decades, nature has been keeping us completely in the dark about the nature of dark matter. 
In view of this, the scope of theoretical and phenomenological studies has become significantly 
broader than it used to be some years ago. An often discussed possibility nowadays is a dark 
sector featuring new gauge interactions. Phenomenologically, such new gauge interactions can for 
instance lead to Sommerfeld enhancement of the dark matter (DM) annihilation cross section Il ¬ 
ls], to DM self-interactions [H [5], which may affect small scale structures in the Universe such as 
dwarf galaxies, or to DM bound states [SHE]. In this context, it is worth noting that there are 
long-standing discrepancies between theory and small scale structure observations, which could be 
explained either by dark matter self interactions BU or by baryonic feedback PEES]. Also recent 
observations of the galaxy cluster Abell 3827 could point towards self-interacting DM El (see, 
however, ED- 

Here, we discuss the prospects of detecting new gauge interactions in the dark sector at the 
LHC. We work in a simplified model with a dark matter fermion x and a U(l)' gauge boson A' 
(dark photon), kinetically mixed with the photon [3j, 13] (see e.g. [H] for a review). This model, 
which we dub radiating dark matter, serves here as a proxy for a large class of more complicated 
scenarios, to which our conclusions can be applied as well. If the dark matter and dark photon 
masses are relatively small compared to the typical partonic center of mass energy at the LHC, 
emission of dark photons from a dark matter particle receives a strong enhancement in the collinear 
direction. This is analogous to the emission of collinear photons in QED or the emergence of a 
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parton shower in QCD. Thanks to the collinear enhancement, DM production at the LHC is 
typically accompanied by the emission of several dark photons whose subsequent decay leads to a 
collimated jet of SM particles. These SM particles will often include electrons or muons, resulting 
in a lepton jet signature. Depending on the dark photon decay length, the SM particles in the 
lepton jets can be produced close to the primary interaction vertex (prompt lepton jets) or further 
away from it (displaced lepton jets). In this paper, we explore both prompt and displaced lepton 
jets from the theoretical and from the phenomenological side. In particular, we derive semi-analytic 
expressions for the spectrum of radiated A' particles, and we compute limits on U{\)' models from 
two ATLAS lepton jet searches mm- We also discuss how the discovery potential will improve 
the 13 TeV LHC. 

Lepton jets have been discussed previously in various contexts. In supersymmetric models, they 
can be a signal of cascade decays of heavy new particles into much lighter states which subsequently 
decay to lepton pairs [EM- Lepton jets could also arise in models featuring non-standard Higgs 
decays to dark sector particles that subsequently decay to leptons [231426] . DM production in 
association with a single radiated A' boson have been discussed in If the A' is heavy, it can 

be reconstructed as a dilepton resonance, while for light A', a narrow hadronic jet from its decay to 
quarks can be a promising signature, especially since an A'-induced jet could be distinguished from 
a QCD jet by using jet substructure techniques [28]. In related works, multi-lepton final states at 
B factories have been studied in [ESI 00]. Several recent papers discuss how strong dynamics in the 
dark sector can lead to the formation of a dark parton shower, with subsequent decay of dark sector 
particles into hadronic jets with non-standard properties like displaced vertices, high multiplicities 
of vertices in the jet cone [31] and significant missing energy aligned with a jet [32]. It is even 
possible that the DM itself interacts relatively strongly with SM particles, so that its production at 
a collider would lead to jets rather than the usual missing energy signature [53]. An experimental 
search for displaced hadronic jets has been presented in [33] by the ATLAS collaboration. 

The plan of this paper is as follows: in sec. [IIJ we introduce our simplified model of radiating 
dark matter and discuss its main phenomenological features. We then analyze the physics of dark 
photon showers in sec. H both analytically and numerically. We derive the expected spectrum of 
radiated dark photons in three different ways: using recursive integral expressions, using an inverse 
Mellin transform to compute the spectrum from its moments, and using a Monte Carlo simulation 
in Pythia 8.2 [351 - 13?] . Our main phenomenological results are presented in sec. IV There, we recast 
the existing ATLAS searches for prompt lepton jets [T5] 


and for displaced lepton jets m to set 
limits on our dark radiation model, and we also compute the sensitivity of future LHC searches at 
13 TeV center of mass energy. We summarize our results and conclude in sec. [V] 


II. MODEL 
A. The Lagrangian 

The results of this paper can be applied to any model in which a < 10 GeV DM particle 
interacts with a new < 10 GeV gauge boson, as long as the typical energy at which DM particles 
are produced at the LHC is much higher than their masses. To illustrate our main points, we 
will use a specific toy model of radiating DM, in which the dark sector consists of a fermionic 
DM particle x and a massive U(l)' gauge boson A'. The coupling strength between \ and A' 
is described by a dark fine structure constant a J = g^,/(47r). Moreover, the dark photon A' is 
kinetically mixed with the SM photon. The Lagrangian for the dark sector is thus 

£dark = X(i0 ~ m x + ig A >A')x ~ ^ F ', l u F ' ,u ' ~ ^m^A^A'" - ^F^F 111 ', (1) 
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FIG. 1: Feynman diagram for dark matter pair production at the LHC through an s-cliannel Z' resonance, 
followed by “dark radiation”, i.e. emission of several— mostly soft or collinear—dark photons. 


where rriA' and m x are the masses of A' and y, respectively. We denote by F^ u and F'^ the 
held strength tensors of the SM photon and the dark photon A' , respectively. The dimensionless 
coupling constant e describes the strength of the kinetic mixing between the A 1 and the photon. 
Due to SU( 2) invariance, the A' should also mix kinetically with the Z. but effects of this mixing 
are suppressed by factors of m A , /M| and are therefore neglected in the following. 

Since the kinetic mixing between A' and the photon is too small to lead to significant DM 
production at the LHC, we assume an additional DM-SM coupling. For definiteness, we take this 
coupling to be through a heavy s-channel vector resonance Z l with couplings to all quark flavors 
and to dark matter. Since the dynamics of dark radiation does not depend on the primary DM 
production mechanism but only on the production cross section and to some extent on the DM 
energy spectrum, our results will apply to any model in which DM particles can be produced in 
significant numbers at the LHC, including for instance models with contact interactions, f-channel 
mediators, or Higgs portal interactions. The relevant terms in the Lagrangian of our toy model are 

Lz' = g q ^2 + 9x x$'x , ( 2 ) 

/ 


where qj is the SM quark field of flavor / and g q , g x are the Z' couplings to quarks and DM particles, 
respectively. For simplicity, we have assumed the coupling to quarks to be flavor universal. 

Fig. a illustrates DM pair production at the LHC via s-channel Z' exchange, followed by radi¬ 
ation of several dark photons. Due to the assumed lightness of y and A ’, the emission is enhanced 
in the collinear direction. Due to this enhancement, with a moderate dark fine structure constant 
a a' ~ 0(0.1), we typically expect a few dark photons to be radiated in each DM pair produc¬ 
tion process. We will discuss the dynamics of this “dark radiation” 
analogous to parton showers in QCD. 


m sec. 


Ill using a formalism 


B. Benchmark points 


We define two benchmark points in the parameter space of our toy model, which we will use 

we will also discuss in detail how 


to illustrate our main points in secs. Ill and IV In sec. IV 


departing from these benchmark points affects our results. The two benchmark points A and B 
are summarized in table [IJ together with several phenomenological observables derived from them. 

In both cases, we assume a Z' mass of 1 TeV. We choose g q and g x such that the resonant Z' 
production cross section is about 1 pb at the 8 TeV LHC for benchmark point A, and about a 
factor of 10 smaller for benchmark point B. In both cases, the Z' has a branching ratio ~ 85% 
for the decay Z l -» yy. We choose both m x and m,A' to be of order GeV or below. Together 
with a moderately large dark fine structure constant otA' = 0.2, this leads to the radiation of 
a significant number of A' bosons when DM is produced at the LHC. Note that the number of 
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[TeV] 



[GeV] [GeV] 


[mm] 

[io- 6 ] 

[Pb] 

[ P b] 

[Pb] 

[GeV] 




A 

1 

0.1 

1 

4 

1.5 

0.2 

10 

2.8 

0.58 

0.85 

2.7 

31.3 

84.8% 

3.50 

3.51 

B 

1 

0.03 0.3 

0.4 

0.4 

0.2 

1 

24 

0.052 

0.076 

0.244 

2.82 

84.8% 

5.15 

5.17 


TABLE I: Values of the model parameters mz> (heavy mediator mass), g q , g x (heavy mediator couplings), 
m x (DM mass), (dark photon mass), a a 1 (dark fine structure constant) and cr (dark photon decay 
length) at our two benchmark points A and B. We also show the resulting values for several derived quantities, 
in particular the kinetic mixing e corresponding to the given ct and m A ', the resonance cross sections a 7 (Z'), 
as(Z') and a 13 (Z') for Z' production pp —> Z' at the 7 TeV, 8 TeV and 13 TeV LHC, respectively, the total 
decay width Tz r of the heavy mediator, its branching ratio to DM pairs, and the average numbers 2 {riA')g 
(2 (?r J 4 ') 13 ) of radiated dark photons in each DM pair production event at the 8 TeV (13 TeV) LHC. 


radiated A' bosons is almost independent of the collider center of mass energy because the boost 
of the primary DM particles, which determines the amount of radiation, is mostly dictated by the 
s-channel Z' resonance. We choose the dark photon decay length cr such that at benchmark point 
A, most dark photons will decay far from their production vertex, while at benchmark point B, the 
displacement of the decay vertex is much smaller. Thus, benchmark point A is in the realm of lepton 
jet searches requiring displaced decay vertices, while benchmark point B is most easily detectable 
in searches for prompt lepton jets. Note that in table |T] we treat cr as a fundamental parameter 
and e as a derived quantity. This will only be relevant when we vary the model parameters around 
their benchmark values in sec. IV (figs. [6] and [7]). To make the discussion there more transparent, 
we choose to keep cr rather than e fixed when varying the other model parameters. Of course, as 
long as all other model parameters are fixed, there is a one-to-one correspondence between cr and 


e. 


Our parameter choices are consistent with existing collider constraints: In particular, dijet 
searches do not exclude Z' bosons with the parameters chosen here [38l 3!) . For instance, the 
ATLAS dijet search in 20.3 fb _1 of 8 TeV LHC data [38] restricts the product a(Z') x BR (Z' —> 
qq) x A of the Z' production cross section, the branching ratio to a quark-antiquark final state, 
and the experimental acceptance factor A ~ 0.6 to be < 1 pb at 95% confidence level (CL). The 
corresponding CMS search |3H] does not explicitly quote a limit at mz' = 1 TeV, but noting that no 
spectral feature is observed at this mass, the limits given on heavier Z' bosons can be extrapolated 
down, indicating that the CMS limit is stronger than the ATLAS limits by perhaps a factor of two. 
The acceptance factor is similar in ATLAS and CMS. Limits on dijet resonances are also provided by 
the CDF experiment at the Tevatron, yielding a(Z') x BR(Z' —)> qq) x A < 0.11 pb [®]. Since CDF 
requires both jets to have a rapidity |y| < 1, we conservatively estimate the acceptance in CDF is 
not larger than in CMS. Including an acceptance factor of 0.6, we predict at our benchmark point A 
(j(Z') x BR(Z' — > qq) x A = 0.08 pb at the 8 TeV LHC and a(Z') x BR(Z' — > qq) x A = 0.001 pb 
at the Tevatron, which is clearly consistent with the above constraints. The pp —>• Z’ —>• XX 
production cross sections at benchmark point B are even smaller. 

DM production through heavy resonances is also constrained by monojet searches mm- 
ATLAS sets a lower limit mz' / y/Sqdx > 2 TeV [JTj, whereas we have mz’ / y/9q9\ ~ 3.2 TeV at 
benchmark point A and mz’ / y/9qSx ~ 10 TeV at benchmark point B. Both benchmark points thus 
satisfy the limit. 

We also need to consider cosmological constraints on the dark sector of our toy model. The 
most important question here is whether x can be a thermal relic and account for all of the DM in 
the Universe. If there is no particle-antiparticle asymmetry in the dark sector, the answer is no. 
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The reason is that at our benchmark points the annihilation cross section |43] 


(<■ TV )xx-* A 'A' ~ 


ira 2 A , (1 -m^/m 2 ) 3 / 2 
m\ (1 - \m? A ,/mty 2 


( 3 ) 


is much larger than the thermal relic cross section of fewx 1CU 26 cm 3 /sec. Reducing the annihilation 
cross section would either require a significantly smaller ct^/, precluding significant dark radiation 
at the LHC, or a larger A' mass, in particular m,A' > m x . th e latte 1 ' case, only the annihilation 
channel yy QQ would remain open. The cross section for this channel is |44l Rf5lj 


/ \ _ „. 3N f9 2 q 9 2 x + m2 q Lrf 

xx^-qq 2ir (4m 2 —fn|,) 2 y m 2 ’ 


(4) 


where Nf is the number of kinematically accessible quark flavors. Since ( av ) X x^qq scales as 
g^g^m^/qn^,, it is much smaller than the thermal relic value for the small x masses that we 
are interested in here. It would thus lead to overclosure of the Universe. For symmetric DM mod¬ 
els, we therefore conclude that y must be a subdominant component of the DM in the Universe 1 . 
This, of course, does not change the LHC phenomenology that we are interested in here, it only 
means that our toy model cannot provide a complete description of the dark sector. The fact that 
for symmetric DM models a discovery of x through a dark radiation signature is only possible if x 
is a subdominant component of DM is in line with the general statement that, in models with mul¬ 
tiple thermally produced DM components, it is likely that the subdominant DM components are 
discovered at the LHC first due to their larger interaction strengths. For asymmetric DM, i.e. DM 
with a primordial particle-antiparticle asymmetry generated at some high scale, x could account 
for all of the DM in the Universe. In this case, the large annihilation cross section to A' pairs is not 
a problem but a feature because it allows for efficient annihilation of the symmetric component of 
the primordial DM soup. However another opposite charge particle might be needed due to gauge 
invariance for such light A' [471j . Let us emphasize that, in this paper, we will indiscriminately 
call x the dark matter, even though the above discussion shows that it does not need to be the 
dominant DM component. 

Next, let us comment on constraints from direct and indirect DM searches. Without a primor¬ 
dial y-y asymmetry, both types of constraints are easily satisfied due to the required small relic 
abundance of y. For asymmetric DM, indirect searches are also insensitive because yy annihilation 
is impossible in such models. The sensitivity of direct detection experiments is limited due to the 
small m x at our benchmark points. For benchmark point B with m x = 0.4 GeV, the DM mass is 
obviously below the direct detection threshold, while for benchmark point A with m x = 4 GeV, 
one might worry about limits from low-threshold experiments. The best limit in this mass range 
comes from CDMSlite and requires the spin-independent y-nucleon scattering cross section <t x n 
to be below 1.5 x 10~ 40 cm 2 if y accounts for all of the DM in the Universe |48j . Our prediction 
at benchmark point A is a x n ~ 6.5 X 10 -42 cm 2 , avoiding the bound by more than an order of 
magnitude. 

A different set of bounds on our model comes from the fact that y particles have self-interactions, 
mediated by the A! boson. Such DM self-interactions are constrained by observations of colliding 
galaxy cluster (most notably the Bullet Cluster) |42j and by measurements of the ellipticity of DM 
halos of groups of galaxies m- The resulting limit on the DM self-scattering cross-section is 

a xx/ m X ~ 1 cm2 /g = 1-78 X 1CT 24 cm 2 /GeV. (5) 


1 There is an exception when A' mass is slightly higher than y mass, one can exponentially reduce the thermal 
averaged annihilation cross section and achieve the right relic abundance 146] . 
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FIG. 2: Branching ratios for the 19 dark photon decay channels included in our analysis as a function of 
the dark photon mass. 


Obviously, this constraint applies only if x forms the bulk of the DM in the Universe. As explained 
above, this is only possible in our scenario if there is a particle-antiparticle asymmetry in the dark 
sector. In the perturbative regime, where aA'm x /mA' < 1, a simple calculation shows that the 
predicted non-relativistic self-scattering cross section is [5] 


a 


/rn x ~ 5 x 10 31 cm 2 /GeV x 


xx 


m. 


Qt4' \ f 

0.1 ) 1 GeV 


GeY V 
m A ' ) 


( 6 ) 


a 


xx 


/ m * 


10 30 cm 2 /GeV 


For our benchmark point A (B) from table [jj this leads to 
(10~ 29 cm 2 /GeV), well within the observational limit. 

Regarding the dark photon mass tua' and its kinetic mixing e with the SM photon, strong 
constraints exist from low energy experiments and from astrophysics. These constraints require 
e < 10 _1 ° for rriA' < 10 MeV, but relax to e < 1CU 3 for 10 MeV < my < 10 GeV, the mass range 
that we are interested in here. Our benchmark points from table |T] are thus consistent with limits 


from dark photon searches. We will discuss these limits in more detail in sec. IV see in particular 

fig. no 


C. Dark Photon decay 

After dark photons have been radiated from pair-produced % particles at the LHC, they decay 
to SM particles via their kinetic mixing with the photon (see eq. 0)- For mA' < 2 GeV and 
via' < 2m x , the dominant decay modes include A! —» e + e _ , / r + ^ _ , as well as decays to various 
hadronic resonances. The partial decay width to leptons of a specific flavor t is 

r(A' - (*f) = UfW ,/l-44 (l + 2-4) . (7) 

3 y mr A , \ mr A , J 

The hadronic decay width can be calculated using measurements of hadron production cross sec¬ 
tions at e + e“ colliders, in particular the ratio R(s) = c(e + e - —> hadrons)/(j(e + e^ —> /r + /Lt _ ), 
where s is the center-of-mass energy. At s = rri A ,. the ratio between the partial A' decay width 
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into a given hadronic final state and the partial decay with to the /r + /x final state should be equal 
to R(s ): 


r(A/ —> hadrons) = r(^4 / -» ^ + /r ) R(s = m\,). 


( 8 ) 


We use hadronic cross-sections from refs. mm to calculate the A' decay branching ratios. We 
include 19 decay channels, shown in fig. [2] In the calculation, care must be taken to avoid double 
counting of hadronic degrees of freedom. For example, the decay A 1 —> 7r + 7r _ 7r 0 7r° receives con¬ 
tributions both from the direct production of pions and from the decay A' —> u;7r°, followed by 
a; —> 7T + 7r _ 7r°. We choose to preferentially include channels with pions and kaons. For heavier 
mesons, we only include decay channels in which there is no double counting. In principle, there 
is also a double counting problem between kaon and pion final states because of the fully hadronic 
kaon decay channels. Fortunately, these channels have very small branching ratios and can be 
neglected. 

For tua' <1 2 GeV, hadronic A' decays are more conveniently described in QCD language as 
decays to quark pairs. The partial decay width to a quark pair of a given flavor qj is 

F(A' -A q f q - f ) = N c Q\ f V(A' -A £+£-) \ mt =m qf , (9) 


where Q qf is the electric charge of qf and N c = 3 is a color factor. At intermediate A' masses 
~ 2 GeV, a very large number of hadronic resonances appears, so that neither the QCD description 
in terms of quark final states nor the low energy effective theory involving just a few hadronic states 
are applicable. We therefore do not consider this mass range in the following, but we will see in 


sec. 


IV that our results on the LHC sensitivity can be smoothly interpolated between the QCD 


regime and the low-energy regime. 


III. DARK PARTON SHOWER 

Before presenting our numerical results, let us provide some analytic discussion of dark photon 
radiation. Part of the following discussion is based on refs. [531154] . We first note that factorization 
theorems ensure that the cross section can be factorized into short range (hard process) and long 
range (shower) contributions [35]. Thus, we can treat the dark photon shower as independent of 
the primary \X production process. The relevant processes in the development of an A! shower are 
X — > x + A' and A' -A- XX- I n the following, we include only the first process, which is shown also in 
fi g- HI Its probability is much larger than the probability for A! —> XX due to the structure of the 
respective splitting kernels. In particular, the splitting kernel for A 1 -A XX is not divergent in the 
soft limit, in analogy to photon splitting in QED. The radiation of gauge bosons from fermions, on 
the other hand, is formally divergent in the soft and collinear phase space regions. In our scenario, 
the dark matter and dark photon masses regularize the radiation probability in these regions by 
providing a natural infrared cut-off. From a phenomenological point of view, we are more interested 
in collinear emission than in soft emission because the decay products of very soft dark photons 
would not be detectable at the LHC. We work under the assumption of strongly ordered emission, 
i.e. the virtuality of the incoming particle in a splitting process is assumed to be much larger than 
the virtuality of the outgoing particles. The latter are therefore treated as on-shell. Within the 
strongly ordered emission assumption, we can consider each splitting x —> X + 7b as an isolated 
process, since a secondary splitting would be a small perturbation with respect to the first one. 

The differential probability for a collinear splitting can be written as 

OLA' . dt . , 

—dxjP x -+ x ( x ,t), ( 10 ) 
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FIG. 3: Radiation of a single dark photon A' from a DM particle %. 


where t is the virtuality of the incoming DM particle and x is the fraction of its energy that 
is transferred to the outgoing DM particle. The quantity P x -> x (x,t) is the splitting kernel that 
encodes the model-dependent details of the splitting probability. We call attention to the fact 
that we treat otA' as independent of t, i.e. we neglect renormalization group running. Since the 
running is logarithmic, it will not play a major role in the development of the dark photon shower. 
Moreover, the running would depend on the full particle content of the dark sector and is thus 


highly model dependent. Eq. (10) can be understood as the squared matrix element of the splitting 


process, multiplied by the propagator of the incoming DM particle and integrated over the phase 
space of the outgoing particles. In the limit t —> 0 (highly collinear splitting), the propagator factor 


1 jt in eq. (10) diverges as expected. It is prevented from going all the way to zero thanks to the 
non-zero dark photon mass. 


A. Kinematics and notation 

We write the off-shell 4-momentum of the incoming DM particle in a given splitting process as 
(see ref. [3] for a slightly different approach) 


P x ,in = (E,0,0,p) . (11) 

Notice that in general E ^ V~§/2. For instance, if the second x particle, which we call the spectator, 
is on-shell and does not radiate, energy momentum conservation requires its three-momentum to 
be (0,0, —p). This fixes the energy of the spectator to be E 2 = p 2 + m x , and hence E > E s + tua' 
if at least one dark photon is radiated in the event. In the following, we will nevertheless assume 
E = E s = y/I/ 2, the so called collinear approximation. Naturally, we expect it to fail when 
tua' + m x ~ or when the opening angle between the splitting products is large, so the splitting 
is not collinear any more. 

The on-shell 4-momenta of the outgoing DM particle and dark photon after the splitting are 
P x ,out = (xE, -kt , 0, yjx 2 E 2 - - m|), (12) 


and 


k = ((1 - x)E, k t , 0, -x) 2 E 2 -k 2 -m 2 A ,) , 


(13) 


respectively. Since the energies of the outgoing particles must be larger than their masses and 
because energies and momenta have to be positive, the following constraints must be fulfilled: 


m x /E < x < 1 — rriA'/E, x 2 E 2 — m x > kf , (1 — x) 2 E 2 — m\, > k 2 . (14) 

The first relation defines in particular the allowed range of the variable x in the splitting probability 
eq. (10). Unfortunately, this range depends on E and thus on all preceding splitting processes. 
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This would preclude us from treating the dark photon shower analytically. Therefore, we will in 
the following take the minimum and maximum values of x to be 


•7-min — m x /E 0 , 2-max — 1 ITlA' / Eq , 


(15) 


where Eq is the initial energy of the DM particle x before emission of the first dark photon. Since 
most emitted dark photons are soft, we expect E to be close to Eq for most splittings. The relations 


(15) would even be exact for m x = = 0, but since we have m x , tua 1 Eq, we expect them to 


be very good approximations also in our case. 
The virtuality t is defined as 


t — ijPx, out T k) 77t x — ^A 1 T 2p Xi out ■ k , 


and the splitting kernel is given by 


P 


x-^-x 


(x,t) = 


1 + x 2 2 (m 2 +m 2 A ,) 


1 — x 


(16) 


(17) 


We see that P x ^ x {x,t) becomes large for x — > 1, where the emitted dark photon is very soft. 
Since we are dealing with massive dark photons, x is kinematically required to always be < 1 (see 
eqs. (14) and (|15[)). In the following, we will neglect the mass dependent term in P x ~^ x and call 
the thus approximated splitting kernel P x ^ x (x). Like the approximation used for x m i n and x max 
above, also this approximation is necessary to make different splittings completely independent of 
each other, a precondition for our analytic calculations. 

In order to also regularize the collinear t 
when choosing the integration limits for t 
that is 


0 divergence in eq. (10), care must be taken also 
The lower bound i m i n for t is obtained when kt —> 0, 


^min(^) — ^A’ 4" 2p x ,out ' 

= 771 ( 4 / + 2^E 2 x(l — x) — \Jx 2 E 2 — m | \J(1 — x) 2 E 2 — m 2 A ,^ , (18) 

while the upper bound is obtained when kt is maximal, within the constraints (14), that is, 

^max(s-) — 277^/ T 2p XiOU t • k , (19) 

with 

^ 7 ,max (x) = min {(1 - x) 2 E 2 - m\, , x 2 E 2 - m 2 x } . (20) 


Interestingly, when x takes its maximum or minimum value allowed by eq. (14), /c^ max becomes 
0, which leads f m ; n = t max and the splitting probability goes to 0. As in eq. (15), we will in 
the following take E —> Eq also when evaluating t m i n and t max . This approximation is necessary 
to make consecutive splitting processes completely independent and thus treat the dark photon 
shower analytically. 


B. Number of emitted dark photons 


A DM particle is produced at the LHC with a certain maximal virtuality f max and then radiates 
in general several dark photons. In each splitting, its virtuality is reduced, until it finally reaches 
the infrared cutoff t m j n . The expectation value for the number of radiated dark photons is 


(n A >) ^ 


Q-A 1 
2t r 


Z ^max /*^max 

dx / 

min " ^min 


j p x^x( x )- 


( 21 ) 
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Note that this compact expression is valid only if we neglect the t-dependence of the splitting kernel 
and of the integration boundaries x' m j n . x max , i m ; n , f max defined in the previous section. 

The probability that a DM particle x radiates exactly m dark photons is given by a Poisson 
distribution ESSh 


Pm = 


~( n A') (n A , 


ml 


( 22 ) 


To see this, note that eq. (10) implies that the probability for no splitting to occur between t max 
and f m i n 


is given by the Sudakov factor 

A(t n 


f 1 = p~( n A') 
i? ‘"max/ — c > 


(23) 


with (riA') defined in eq. (21). The probability for exactly one splitting is then 


Pi = 


ClA' 

2vr 


<. r imax 

dx / ~^(^miri5 ^max) 


= e ~( n A') 


{nA') ■ 


(24) 


In the first line, the factors A(f m ; n ,t) and A(f,f max ) give the probabilities for no splittings to 
happen in the intervals [t m i n ,t) and (t,t max ], respectively. The splitting kernel P x -^ x {x) describes 
the probability that a splitting happens at virtuality t. For two splittings, we have in analogy 


P 2 = 


A' 

2vr 


dx 


1 ( OLA' 
ll 

~( n ai) M 2 


~ p-( n A')-Z- 

2 ! V 2vr 


* dt 
t 

dx 


^min 
C tmai 


c ft dt r 

dx' / —A(t min ,t)P x ^ x (x)A(t,t')P x ^ x (x')A(t',t ma , x ) 
Ju ■ 1 


dt 


t/ min 

1-max 


c r^max fj-f-t 

dx' / ^-P x ^ x {x)P x ^ x (x') 


= e -( n A') 


2 ! 


(25) 


In the second line, we have extended the integration region of the t' integral to the full range 
[Imin- imax] and included a factor 1/2! to compensate for the double counting induced that way. 
This is the place where the approximations discussed in sec. Ill A are crucial: extending the 
integration interval is only possible if we can assume that x m ; n , x max , f m ; n and t max are the same 
for all splittings and that P x ^ x is independent of t. Eq. (25) straightforwardly generalizes to larger 
numbers of splittings, yielding the Poisson probabilities p m in eq. (22). 


C. Recursion relation for the \ and A' energy spectra 


Let us now study the energy distributions of x anc l ^ in events with radiation of multiple A' 
bosons. Our starting point is the observation that, in the collinear limit, each splitting process is 


independent. As explained in sec. Ill A, we will make the crucial approximation that the integration 
boundaries x m i n , x max , i m in, tmax are the same for each emission process and are given by eqs. (15), 
(18) and ( |19| ) (with E -A Eq). We moreover neglect again the f-dependence of the splitting kernel 
P x -+ x - These approximations will allow us to compute the energy distributions recursively. 

If we call the final energy of the on-shell DM particle after final state radiation E x and the 
energy fraction retained by the DM particle after all splittings X = E x /E q, 2 we can write the final 


2 We use capital letters to denote quantities normalized to the initial energy Eq of the DM particle before the first 
splitting, and lower case letters for quantities normalized to the energy of the DM particle as it enters a particular 
splitting process during the development of the A' shower. 
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DM energy distribution as 


f x (X)=J2Pmfx,™(X), (26) 

m =0 


where f X)7n (X) is the energy distribution of x i n events with exactly m emitted dark photons. 
f Xt m{X) obeys the recursion relation (see also [57] ) 

f l V\ f maX J™ f \ fx, m (X/x m ) s' V s' „ \ ( c \ r 7\ 

Jx, m+l (-^0 — / Jx, l\Xm) ^ X — ^max) ? (27) 

** ^min ^TYl 

where 

1 CV M /*^ m ax 

/xdPO = w—y— / —P x ^ x (X) 0(x m in < X < X max ) , (28) 

is the energy fraction for dark radiation showers with only a single splitting and 0( •) denotes a 
window function that is 1 when the condition in parentheses is satisfied and 0 otherwise. Note that 
for m = 1, the definitions of x as the fraction of energy the DM particle retains in a single splitting 
and of X as the fraction of energy the DM particle retains after all splittings are identical. For 
m = 0, we define f Xj o(X) = 5(1 — X). Note that f XtTn {X) is normalized according to the condition 


JdXf x , m (X) = 1. (29) 

When evaluating f x (X) numerically, we truncate the recursive series at m = 10, which is justified 
by the fact that the Poisson probability p m is very small at m ^$> (n A i). 

The distribution of the energy E A t of radiated A! bosons can be obtained in a similar way. With 
the notation Z = E A i/Eq, it is given by 

- oo m 

f A '{z) = 7-^ E > (30) 

' A '' m =1 k= 1 

where f A \k( z ) is the energy distribution of the k- th emitted A' boson. Note that f A \k( z ) is 
not a function of the total number of radiated dark photons m because each splitting process is 
independent. The recursion relation for f A 'k{Z) is 

fA',k+i ( z ) = [ dx k f X) l(.Tfc)^ 4 ' k ( Z / Xk ^ 0(1 - x' max < Z < 1 - x min ). (31) 


D. Energy spectra from an inverse Mellin transform 

A particularly elegant way of computing the energy spectrum f x (X) of dark matter particles 
and the spectrum f A >(Z) of dark photons is by first computing the moments (X s ) and ( Z s ) of 
these distributions, and then applying an inverse Mellin transform to recover f A >(Z) itself. The 
Mellin transform of a function f(X) is defined as [5Sj 

n oo 

M[f](s) = p(s) = dXX s ~ 1 f(X), (32) 

Jo 

and the inverse transform is given by 

i l>c+ioo 

f(X) = — dsX- s <p(s). (33) 

J c —200 
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; dt 


The Mellin transform and its inverse are closely related to the Fourier and Laplace transforms [58] 
and can therefore be evaluated efficiently using the FFT (Fast Fourier Transform) algorithm oT)h 
For events in which a single A' is emitted, the moments of the DM energy distribution after the 
emission, weighted by the probability of having exactly one emission, are given by 

r\j At f ^ ma x r ^max 

Pl{X') 1A ,=e-M^ dxx' / 

= e{n A ')X*. 

In the last line 

a a 1 /'' Tinax , s f tmax dt 




‘■min 


\-^x( x ) 


we have defined 

x 1 1 


(34) 

(35) 


r X max ft 

7 - r — / dxx s / —P*_~(x). 

{n A i) 2vr J x . J t . t x ^ xK J 

As in the previous section, we again take x' m ; n , x max , i m i n and t max to have the same value for 
each splitting process, and we neglect the f-dependence of Px-^x- If two A' bosons are emitted, we 
obtain analogously 




^max 


dx x s 




'min 

'^max 


rt max max 

l t . T 

‘-min 


ft (Jft 

dx'x' s / v P x ^ x (x)P x ^ x (x f ) 

” ^min 


®min ” ^min ” ‘■min 

, /' tmax dt f Xinax d / is f tmax dt' 


- e ~( n A') ( n ' A '} X s2 

2! 


dx x s , 

Jt ■ t 

min ‘min 


"‘'min 


/“‘max 

I T P x ^ x {x)P x ^ x {x') 

^min 

(36) 


In the second line, we have again used the approximations from sec. Ill A to extend the integration 


region of the t' integral (see the discussion below eq. (25)). Eq. (36) can be easily generalized to 
the case of m emitted A' bosons: 


Pm (X S 


>mA ml 


Summing over m, we find the moments of the overall A! energy distribution: 


ip(s + l) = (X s ) = J2p m (X s 


' mA' 


= e ~(n A i)( 1-X‘) 


(37) 


(38) 


m =0 


An inverse Mellin transform of <p(s) then yields the dark matter energy spectrum f x (X). As can 


be seen from eq. (33), this requires evaluating y>(s) at complex values of s. A numerically stable 


way of doing this is by using the FFT algorithm to compute X s (eq. (35)), then using eq. (38), 


and afterwards apply another FFT to evaluate the inverse Mellin transform. 

The procedure for obtaining the A' spectrum is alike and yields for the contribution of showers 
with m dark photons to the moment s, weighted by the probability of emitting exactly rn dark 
photons, 


Pm {Z s 


mA' 


(n A i) ml 


with 


v? - 1 

( n A >) 2 tt , 

Summing over all emissions we obtain 

( z s ) = 


dx (1 — x) s 


Z s 1 — e -( n A')P~ Xs ) 


-z*Y, xak ~\ 

k=l 

(39) 

■tmax ,U 

t l J x- x >x( x ) • 

(40) 


(n A i) 1 - 

An inverse Mellin transform yields the A! spectrum f A ’{Z). 


(41) 
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m x 

[GeV] 

m A ' 

[GeV] 

2(n A >) 

2 (riA' ) Pyt hia. 

(X) 

(X) Pythia 

(Z) 

(Z) Pythia 


50 

1.5 

2.130 

2.340 

0.873 

0.837 

0.119 

0.140 


50 

0.4 

2.848 

3.084 

0.871 

0.835 

0.091 

0.107 

A 

4 

1.5 

3.476 

3.540 

0.729 

0.697 

0.156 

0.171 


4 

0.4 

4.990 

4.825 

0.712 

0.681 

0.116 

0.132 

B 

0.4 

0.4 

5.691 

5.215 

0.626 

0.608 

0.132 

0.150 


TABLE II: Characteristics of the dark photon shower for various choices of the DM mass m x and the dark 
photon mass m^. The rows labeled “A” and “B” correspond to the two benchmark points from table [T} In 
all cases, we assume x pah production at a center of mass energy y/I = 1 TeV and we take ola' = 0.2. We 
show the predicted number 2 (ua 1 ) of dark photons per event (with the factor of 2 coming from the fact that 
we consider \ pair production), the average fraction (X) that a \ particle retains of its initial energy after 
showering, and the average fraction (Z) of the initial \ energy that each dark photon receives. We compare 
the predictions of our semi-analytic calculations, eqs. ©, ( [38] ) and pTj ) (with s = 1), to the results of a 
Monte Carlo simulation in Pythia. As expected, the results satisfy the relation ( X ) + ( ua>) (Z) = 1. 


E. Comparison of analytic results to Monte Carlo simulations 


We have also simulated the radiation of A' bosons fully numerically in Pythia 8 [35143*7] . using 
the “hidden valley” model implemented therein. Note that, like our analytic calculations, also the 
dark photon shower implemented in Pythia includes only dark photon radiation, x and 

X —> A'dP, not dark photon splitting, A' —> XX- 

To compare the x and A! energy spectra obtained using the different approaches, we assume 
X pair production at a center of mass energy of \f% = 1 TeV. In Pythia, this is achieved by 
simulating the process e + e _ — >• Z’ — > xXi followed by the dark photon shower. The results for 
several characteristic parameters of the dark parton shower, namely (ua 1 ), (V) and {Z), are shown 
in table [TI} C omparing Pythia results to the results obtained using eqs. (21), eq. ( p8| ) (with s = 1) 
and eq. (41) (with s = 1), we find excellent agreement. We attribute the remaining differences to 


the approximations we had to make in the analytical calculation, in particular the assumption that 
the upper and lower integration limits in x and t do not depend on x and t themselves. 

A more detailed comparison between our analytic results and Pythia is shown in figs. [4] and [5} 
In fig. [4j we have plotted the distribution of the number of emitted dark photons per event (i.e. 
per XX pair) at our two benchmark points from table B Fig. E shows the energy spectra of x and 
A' after the dark photon shower has developed. 

As we can see in fig. [4j Pythia results confirm that dark photon emission is approximately a 
Poisson process. We attribute the small deviations to the approximations we made in deriving 
eq. (22) in sec. Ill B 


Regarding the energy spectra plotted in fig. [5[ we find excellent agreement between both of our 
analytic approaches (red dashed and blue dot-dashed curves) and also with Pythia results (black 
dashed curves). Note that the recursive method is numerically less stable than the Mellin trans¬ 
form approach because numerical errors—especially those incurred close to the upper integration 
boundary x max at which the splitting kernel has a regularized singularity—can accumulate in the 

We have also carried out simple leading order simulations of the 

(gray solid curves in fig. [ 5 ]) 


recursive formulas (27) and (31) 

processes e + e - —> Z' — > XX and e + e _ —> Z' — > xxA' i R CalcHEP 
and find that, as expected, this simulation reproduces our analytic result for the specific process 
e + e~ —> Z l —> xxA' very well, but is insufficient as a proxy for the full dark photon spectrum. 
Since the full dependence of the matrix element on m x and tha' is taken into account in CalcHEP, 
it also illustrates that the approximations discussed in sec. 


Ill A are justified at this level. In 
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FIG. 4: The distribution of the number n of dark photons emitted in each \X Tail' production event at a 
center of mass energy of Vs = 1 TeV. The model parameters are given in table [ij The solid curves labeled 
“analytic” show the Poisson probability e~ 2< ' riA '' > [2 (riA')] n /n\, with (n^) given by eq. (211. The factors 
of 2 arise from the fact that two DM particles are produced in each event. The dotted curves show the 
distribution obtained from a Monte Carlo simulation of the dark photon shower in Pythia. 


fig. [5j we also show the contributions p m fx,m{X ) and ( 1 / (ua 1 )) Ylk=i PmfA',k{Z) from events with 
a particular fixed number of dark photons to the overall energy spectra f x (X) and /q/(Z) (dashed 
rainbow-colored curves). Note that events with several emitted A! bosons can dominate over events 
with a single A' because 2 {tla') > 1 at our benchmark points. This illustrates that the A! spectrum 
would be very difficult to obtain from a full matrix element calculation and necessitates a parton 
shower algorithm. 

In summary, we have developed in this section two analytic methods for computing the energy 
spectrum of x an d A\ and we have verified them by comparing to Monte Carlo simulations. 


IV. COLLIDER SEARCHES FOR LEPTON JETS FROM DARK RADIATION 


7 TeV and 
“hidden valley' 


To constrain radiating DM models using existing LHC lepton jet searches, we numerically 
simulated the expected signal and background event spectra at center of mass energies of ^/s = 
s = 8 TeV. As in sec. Ill E, signal events are generated using the implementation of 

We compare our simulation results to data from a 
15] and from an 8 TeV ATLAS search for displaced 


a "mclclen valley" model in Pythia 8 
7 TeV ATLAS search for prompt lepton jets 
lepton jets m- We have also simulated events at y/s = 13 TeV to predict the sensitivity of future 
searches for prompt and displaced lepton jets. 


A. Prompt lepton jets 

Our prompt lepton jet analysis follows the 7 TeV, 5 fb” 1 ATLAS search in ref. 115) . but including 
only muonic lepton jets. The reason for discarding lepton jets with electrons is that it is very 
difficult to implement the corresponding selection criteria without a full detector simulation. The 
trigger requirement for muonic lepton jets is the presence of either three muons with transverse 
momenta pt > 6 GeV or one muon with px > 18 GeV. Muons are required to have pseudorapidity 
\r]\ < 2.5, and they are required to have a track in the inner detector associated with them. We 
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(a) 


(b) 




(c) 


(d) 


FIG. 5: (a), (b) Energy spectrum f x (X) of DM particles x after final state radiation; (c), (d) energy spectrum 
fA'(Z) of dark photons A' emitted as final state radiation. The panels on the left are for benchmark point 
A from table [TJ the panels on the right are for benchmark point B. In all cases, we assume \ pair production 
at a center of mass energy \f§ = 1 TeV. We compare the results from the recursion formulas described 
in sec. |fflC| the Mellin transform method developed in sec. |TTTD] the dark photon shower simulation in 
Pythia, and a simple leading order simulation of e + e - —>■ Z' —> XXA' in CalcHEP. For the Mellin transform 
method, we also show the result separated according to the number of A' bosons emitted in each XX pair 
production event. 


implement the last requirement by demanding that the parent A' of the muon decays within a 
radius of 122.5 nun (corresponding to the position of the last silicon pixel layer) from the beam 
axis. Finally, muon tracks are required to have a transverse impact parameter |do| < 1 mm with 
respect to the primary vertex. 

We cluster muons into lepton jets using the following algorithm: we collect all muons with a 
cone of radius A R = 0.1 around the highest-p^ muon in the event to form the first lepton jet. We 
then iterate this procedure, starting, in each iteration, with the highest-p^ muon not associated 
with a lepton jet yet. Ref. m then defines two classes of events that are considered in the analysis: 

Double muon jet events require at least two muonic lepton jets, each of which must contain at 
least two muons with px > H GeV. If the event was selected only by the single muon trigger, the 
leading muon in the lepton jet is in addition required to have px > 23 GeV. Moreover, the two 
muons closest in px within each lepton jet are required to have an invariant mass < 2 GeV. 
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To reject non-isolated lepton jets, the variable 


P 


= Ei e ta 
Pt, lj 


(42) 


is defined, where the sum in the numerator runs over the Ex of all calorimeter deposits within 
A R = 0.3 of any muon in the lepton jet, but excluding deposits within A R < 0.05 around any 
muon. The denominator is the px of the lepton jet. The isolation cut imposes p < 0.3. 

Single muon jet events require only one muonic lepton jet, but this lepton jet must contain at 
least four muons with px > 19 GeV, 16 GeV, 14 GeV for the three leading muons and px > 4 GeV 
for all other muons. The same invariant mass cut as for the double muon jet analysis is imposed, 


and the isolation variable (42) has to satisfy p < 0.15. 


Since we find that double muon jet events are much more sensitive to our radiating DM model 
than single muon jet events, we do not include the latter in the following analysis. 

The dominant background to the muonic lepton jet search is from misidentified QCD multijet 
events, which contribute 0.5 ± 0.3 events in the signal region m- For the 7 TeV analysis, we 
can use this number at face value, while for the extrapolation to y/s = 13 TeV, we assume the 
background cross section after cuts to be a factor of 3 larger. This factor of 3 corresponds roughly 
to the scaling of the cross section for QCD jet production m- We assume the relative error on 
the background to remain the same as at 7 TeV. Note that, to be able to use this simple scaling 
of the background cross section, we have to use the same cuts at 7 TeV and at 13 TeV. 


B. Displaced lepton jets 

Our search for displaced lepton jets, based on the ATLAS 8 TeV, 20.3 fb” 1 analysis [T6j, is 
somewhat more involved than the search for prompt lepton jets described in the previous section. 
First, we need to be more careful in simulating detector effects, which we do at particle level. 
We pay particular attention to the displaced A' decay vertices, which lead to unusual detector 
signatures for each particle. The coordinates of the A' decay vertex and the momenta of the A' 
decay products are smeared by Gaussians with widths proportional to 1/L xy , where L xy is the 
transverse distance of the A' decay vertex from the beam pipe. This way, we in particular take 
into account the fact that the momenta of charged particles that travel through more material are 
affected more strongly by the smearing. All parameters of our detector simulation are calibrated 
against the lepton jet reconstruction efficiencies shown in fig. 6 of ref. m • Since this plot is for 
a sample of dark photons with flat px and r] distributions in the range px £ [10,100] GeV and 
r/ E [—2.5, 2.5], we generated a similar event sample for the comparison. To further improve the 
agreement between our predicted efficiencies and the results from ref. m, we apply a fudge factor 
in each L xy bin. 

Following m, muons are selected from the event sample if they are produced inside the sensitive 
ATLAS detector volume, i.e. up to a radius of 7 m. Only loose requirements on the transverse 
momentum and pseudo-rapidity of reconstructed muons are applied: their px must be above the 
trigger threshold of 6 GeV and the pseudorapidity \rj\ must be < 2.5. Since the goal of the search 
in ref. m was to constrain long-lived dark photons with displaced decay vertices, the background 
from prompt muons is reduced by requiring stand-alone muons, i.e. muons that are detected in the 
muon spectrometer but cannot be matched to a track in the inner detector. This essentially means 
that the dark photon decay in which the muon was produced must have happened outside the last 
pixel layer of the inner detector (ID), located at a radius of 122.5 mm from the beam pipe. To 
suppress cosmic ray muons, the reconstructed muon trajectory must still be associated with the 
primary interaction vertex. In particular, it must fulfill requirements on the impact parameters 
\zq\ < 270 mm (longitudinal) and |do| < 200 mm (transverse). 
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Hadronic jets (or, more precisely, calorimeter jets) are based on energy deposits in the elec¬ 
tromagnetic calorimeters (ECAL) and hadronic calorimeters (HCAL), clustered using the anti-fc^ 
algorithm [62] implemented in FASTJET {63] with a cone radius of R = 0.4. Calorimeter jets must 
fulfill pt > 20 GeV and \rj\ < 2.5. 

Three different types of lepton jets (LJ) are defined: A type-0 or “muonic” LJ arises from dark 
photon showers in which all A' bosons decay to muons. It consists of at least two displaced and 
collimated muons with no additional activity nearby. If however, such muons are accompanied by 
a single calorimeter jet, we call the sum of these objects a type-1 or “mixed” LJ. A type-2 LJ or 
“calorimeter” LJ is defined as an isolated calorimeter jet that fulfills additional requirements to 
distinguish it from QCD background. Note that, especially for type-1 and type-2 LJs, the term 
lepton jet used in [T6] is somewhat misleading because these objects can also arise from hadronic 
A' decays. The clustering algorithm for all three LJ types is described in the following. 

We collect all muons and calorimeter jets inside a cone of radius A R = 0.5 around the highest- 
px muon in the event. If at least one other muon and no calorimeter jet is found inside this cone, 
those muons are combined to a muonic LJ. If at least one muon and exactly one calorimeter jet is 
found inside the cone, all objects are combined into a mixed LJ. If no second muon is found within 
A R = 0.5 of the highest-p-p muon, the muon is discarded. The algorithm then continues with the 
highest-pT muon that has not been associated with a lepton jet or discarded yet. After all muons 
have been processed, the remaining calorimeter jets are reconstructed as calorimeter LJs if their 
electromagnetic (EM) fraction is lower than 0.1. This is for example true for dark photons decaying 
to charged pions, but also for a decays to electrons happening inside the hadronic calorimeter. 
Since the pseudo-rapidity region 1.0 < |?y| < 1.4 tends to create a fake low EM fraction due to 
the transition between the barrel and endcap EM calorimeters (ECAL), this region is excluded. 
Additionally the width IE of a calorimeter LJ, defined by 


W = E. Ag -Pr < 0 1, 


(43) 


needs to be small. Here, the sum runs over all particles in the lepton jet. An isolation cut is 
imposed on all reconstructed lepton jets by requiring that the scalar px sum of all charged tracks, 
reconstructed by the inner detector inside a cone of radius 0.5 around the lepton jet axis, is below 
3 GeV. Only charged tracks with px > 400 MeV and with impact parameters |zo| < 10 mm and 
| do I < 10 mm are considered in this procedure. An event is discarded if there are not exactly 
two lepton jets fulfilling \Acj)\ > 1.0. The requirement of two lepton jets leads to the 6 different 
combinations 0-0, 0-1, 0-2, 1-1, 1-2, and 2-2, where e.g. 0-1 corresponds to an event reconstructed 
with one type-0 (muonic) and one type-1 (mixed) lepton jet. 

The sensitivity of the displaced analysis will depend strongly on the A' decay mode. Therefore, 
we list in table III the most important A' decay modes and their associated signatures. Check 
marks (/) indicate in which region of the detector a particular decay has to happen in order to 
potentially be reconstructed as a displaced lepton jet. For decays that would fail our cuts, we 
give the reason why they are vetoed. For instance, due to the requirement of a low EM fraction, 
A! —> e + e _ is only identified as a displaced lepton jet if the decay happens inside the hadronic 
calorimeter, while A' —> p + p~ is also selected if the decay happens inside the electromagnetic 
calorimeter. Decays to any final states involving charged particles are vetoed if they happen in the 
inner detector. Regarding decays to mesons, it is important to note that 7r , K± and K £ have 
lifetimes of order 10~ 8 sec |5U and are thus stable on detector scales (i.e. they will be stopped in 
the hadronic calorimeter before they decay). Neutral pions, on the other hand, decay very fast. 
For Jig, the lifetime of 9 x ICG 11 sec may allow it to travel for several cm before decaying (mostly to 
7r + 7r _ or 7r°7T 0 ). Therefore, A 1 —> tt + tt~ , K + K~ are accepted if they happen in the electromagnetic 
or hadronic calorimeter. The decay A' —> 7r + 7r _ 7r°, on the other hand is vetoed if it happens in 
the electromagnetic calorimeter because the quasi-instantaneous decay of the 7r° to two photons 
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Detector 

A' —>■ e + e _ 

A! ->■ p+p~ 

A! — \ 7T + 7T - / K + K~ 

A' TT + TT~TT 0 

A' -> K° l K% 

LJ type 

2 (calorimeter) 

0 (muonic) 

2 (calorimeter) 

2 (calorimeter) 

2 (calorimeter) 

ID 

track 

track 

track 

track 

U) 

ECAL 

EM fraction 

/ 

/ 

EM fraction 

(/) 

HCAL 

/ 

/ 

/ 

/ 

/ 


TABLE III: Illustration of where in the detector a specific A! decay must happen in order to potentially 
be reconstructed as a lepton jet in our displaced lepton jet analysis (see sec. IV B and ref. |16j. For decays 
that will be vetoed, a reason for the veto is given. The type of lepton jet as which each decay mode is most 
likely to be reconstructed is given at the top of the table. 


deposits a large amount of energy, violating the cut on the EM fraction. The decay A' —> K^Kg is 
successfully reconstructed as a displaced lepton jet in many, but not all, cases. It will be vetoed in 
the inner detector if both the A' and the secondary Kg decay very fast (well within the 12.25 cm 
radius of the inner detector) in spite of their typically large boosts, and the Kg decay mode is 
to 7T + 7r - (branching ratio 69%). It will be vetoed by the cut on the EM fraction if the A' decay 
happens in the inner detector or the electromagnetic calorimeter and the Kg decay mode is to 
7r°7r° (branching ratio 31%). 

is most likely to be 


Note that, as can be read from table III the decay mode A' —> p + p 


reconstructed as a type-0 (muonic) lepton jet, while all other decay modes typically lead to type-2 
(calorimeter) lepton jets. Type-1 (mixed) lepton jets are obtained only if several A' are radiated 
from the same dark matter particle, with at least one of them decaying to p + /r~ and at least one 
of them decaying to a different final state. Therefore, the requirement for a type-1 (mixed) lepton 
jet are the most difficult to satisfy. The predicted signal event numbers in table. IVl confirm this. 


C. Results 


In tables IV and [V] we present the signal and background predictions as well as the observed 
event numbers for the two benchmark models defined in tabled Table IV is for the prompt lepton 
jet analysis from ref. [15] . using 5 fb -1 of 7 TeV ATLAS data (see sec. IVA, while table |v| is fo r 
the displaced lepton jet search from ref. [16], using 20.3 fb -1 of 8 TeV ATLAS data (see sec. IVB). 


Note that, for the prompt search, we include only muonic lepton jets, which are easier to model 
without a full detector simulation. For the displaced search, we show event numbers with and 
without inclusion of the background-limited sample of events with two type-2 lepton jets. We 


see from tables IV and [V] that, as expected, benchmark point A is more easily detectable in the 
displaced search, while benchmark point B is probed more sensitively by the prompt search. 

We also show our predictions for 100 fb -1 of 13 TeV data. In general, we expect better sensitivity 
at 13 TeV because of the larger integrated luminosity expected. Since an estimation of the multijet 
background at 13 TeV is problematic, we restrict the displaced lepton jet search at 13 TeV to type-0 
(muonic) LJs, for which QCD backgrounds do not play a major role. The only relevant source of 
background in the 13 TeV displaced search are then cosmic rays, whose flux is independent of the 
collider center of mass energy as long as we use the same cuts as at y/s = 8 TeV. Our predictions 
for the sensitivity of a displaced lepton jet search at 13 TeV LHC are thus very conservative. For 
the prompt search at 13 TeV, we scale the background cross section from the 7 TeV analysis by a 
factor of 3. 

We use the CLs method [165] to compare our predictions at y/s = 7 TeV and y/s = 8 TeV to the 
data and to constrain the production cross section for a DM pair, a(pp —> Z')BR(Z' xx)- We 
use the ATLAS background predictions p~5l fl6] as the null hypothesis, and the sum of the ATLAS 
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7 TeV 

13 TeV 

Benchmark A 

0.8 

109 

Benchmark B 

3.9 

334 

All backgrounds 

0.5 ±0.3 

30 ± 18 

data 

3 



TABLE IV: Predicted event rates for the prompt lepton jet analysis (see sec. IV A and ref. [15]) at the 


benchmark parameter points from table [T] compared to the background predictions and the observed event 
rates from ref. 1151. 







Lepton jet type 




0-0 

0-1 

0-2 

1-1 

1-2 

2-2 

All 

All excl. 2-2 

Cosmic ray bkg. 

15 

0 

14 

0 

0 

11 

40 ± 11 ±9 

29 ± 9 ± 29 

8 TeV 

Multi-jet bkg. 

Benchmark A 

14 

3 

104 

0 

14 

70 ± 58 ± 11 

200 335 ± 18 ± 100 

12 ± 9 ± 2 

135 ± 12 ±41 

Benchmark B 

2.1 

0.4 

3.0 

0 

0.3 

1.2 

7 ±2.1 ±2.6 

5.8 ±1.7 ±2.4 

data 

11 

0 

11 

4 

3 

90 

119 

29 

13 TeV 

Benchmark A 

Benchmark B 

169 

28 









TABLE V: Predicted event rates for the displaced lepton jet analysis (see sec. IV B and ref. usd at the 
benchmark parameter points from table [T] compared to the background predictions and the observed event 
rates from ref. [16] . In the last two columns, the first error is the statistical uncertainty, while the second 
one is systematic. Our sensitivity study at \J s = 13 TeV includes only type 0-0 events (only muonic lepton 
jets) because a reliable extrapolation of the multijet background to 13 TeV is difficult. 


background and our signal prediction as the test hypothesis. For y/s = 13 TeV, we compare our 
signal+background prediction to the background-only prediction to obtain the expected limit. We 
assume a systematic uncertainty of 30% on the signal normalization to account for errors in the 
signal and detector simulation. The results are presented in figs. [6] and [7] for the two benchmark 
points from table [TJ where in each panel one of the model parameters is varied while the others 
are kept constant. For the displaced search at y/s = 8 TeV, two different limits are calculated, 
one including all tagged events (solid red curves) and one excluding type 2-2 events, i.e. events 
with two type-2 lepton jets (solid black curves). Due to the larger background in the 2-2 sample, 
this second limit is stronger in certain parameter regimes than the limit including all lepton jets. 
For the displaced search at y/s = 13 TeV, (red dotted curves), we use only events with two 
type-0 lepton jets, as explained above. We also show in figs. [6] and [7] the predicted values for 
(j(pp —> Z / )BR(Z / —> xx) as thin dotted curves. 

Let us first discuss fig. pi corresponding to benchmark point A. (We will comment on the 
differences compared to fig. uTbelow.) Since for fixed m^, the dark photon decay length cr depends 
only on the kinetic mixing parameter e, figs, [b] (a) and (b) are equivalent. We see that the best 
limits from the prompt search are obtained at cr < 1 mm (e > 10~ 5 ), while for larger decay length, 
the cut on the impact parameter |do| < 1 mm restricts the sensitivity. For the displaced search, 
the sensitivity is optimal around cr = 10-100 mm, corresponding to e ~ few x 10 -6 . For shorter 
lifetimes (larger e) the dark photon decay vertices tend to be closer to the primary interaction 
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FIG. 6: 95% CL upper limits on the DM DM + nA' production cross section as a function of the model 
parameters for benchmark point A from table [I] In each panel, we vary one parameter while keeping the 
others fixed at their benchmark values. Exclusion limits from the 7 TeV ATLAS search for prompt lepton 
jets jT5| (solid blue) and from the 8 TeV ATLAS search for displaced lepton jets [TB] are shown. For the latter 
search, we show results including all lepton jet events (red solid) and excluding events with two type-2 lepton 
jets (black solid). The predicted sensitivity of similar analyses at V s = 13 TeV is shown as blue/red dotted 
curves. The black dotted lines in each panel show the theoretically predicted production cross sections. 
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FIG. 7: Same as fig. [6j but for benchmark point B from table [T| 
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point, so that the A! decay products are more likely to leave tracks in the inner detector, thus 
failing the displaced lepton jet selection criteria. For too large lifetimes (smaller e), on the other 
hand, most dark photons will decay outside the ATLAS detector and not lead to a signal at all. 
The exact position of the most sensitive regime depends strongly on the A' kinematics and its 
mass. 

For the displaced search, notice also the slight shift of the most sensitive point in cr and e 
between the red solid, black solid and red dotted curves in figs. [6] (a) and (b). This is related 
mostly to the inclusion/exclusion of different types of lepton jets in the three cases. In particular, 
as explained above, we include only type-0 (muonic) lepton jets in the 13 TeV analysis. The decay 
A' —> p + p~ can be successfully reconstructed for shorter A' decay lengths than many other A! 
decay modes (see table |III| ) Consequently, the exclusion limit based on displaced type 0-0 events 
only (red dotted curves in fig. [6] (a) in (b)) is best at relatively small cr (relatively large e), followed 
by the sample excluding type 2-2 events and by the sample including all types of lepton jets. 

Similar arguments help us understand the exclusion limits as a function of the A' mass in fig. [6] 
(c). The different peaks and dips in this plot reflect the dependence of the A' decay branching 
ratios on rriA' illustrated in fig. [ 2 J see also [(3j. The most notable features arise due to the broad p° 
resonance around 780 MeV with the narrow cj resonance on top of it and the narrow 0 resonance at 
1 GeV. The p° resonance decays almost exclusively to 7r + 7r _ pairs, so that around this resonance, 
the relative importance of all other decay channels is diminished. This implies that the prompt 
lepton jet search, which includes only muonic lepton jets, has a greatly reduced sensitivity around 
the p resonance. In the displaced search, most A! —> p° —> 7r + 7r~ decays are reconstructed as 
type 2 (calorimeter) lepton jets (see table III), therefore the sensitivity improves around the p° 
resonance if type 2-2 events are included in the analysis (red solid curve in fig. ©(c))- If type 
2-2 events are excluded (black solid and red dotted curves in fig. [6] (c)), the sensitivity decreases. 
The dominant decay of the u resonance is to 7r + 7r _ 7r°. This again decreases the sensitivity of the 
prompt search around the resonance. For the displaced search, we can read from table III that the 
decay A 1 —> u —> 7r + 7r - 7r° is more likely to be vetoed than A' —> , therefore the sensitivity is 

reduced at the u resonance compared to the off-resonance region. The behavior at the cf) resonance 
is similar to that at the p° resonance. The prompt search is insensitive because the cf) mostly 
decays hadronically. For the displaced search, detection prospects for the dominant decay channel 
4> —> K + K~ (branching ratio 48.9%) as a type-2 (calorimeter) lepton jet are very similar to those 
of the p° resonance decaying to 7r + 7r~ (see table III). The secondary decay 1 —> K®Kg (branching 
ratio 34.2%) has a high probability of being reconstructed as a displaced type-2 lepton jet even if 
the A' decay happens in the inner detector. For these reasons, the sensitivity at the </> resonance is 
further boosted in analyses including displaced type-2 lepton jets. For A' masses above 1 GeV, the 
sensitivities of both the prompt and the displaced searches decrease rapidly. For the prompt search, 
the invariant mass requirement m^ < 2 GeV restricts the realm of sensitivity to uia 1 < 2 GeV. 
For the displaced search, the smaller A 1 boost factor at larger tua 1 implies that the A! decay length 
becomes shorter and events are more likely to leave tracks in the inner detector and be vetoed. 
Moreover, the average number of radiated A' decreases at larger As discussed in sec. [TlJ we do 
not consider the region 1.7 GeV < < 2.3 GeV, where numerous hadronic resonances appear 

but the quark model is not applicable yet. In fig j6] (c), we have smoothly interpolated through 
this region. As can be seen, the low-energy description in terms of hadron final states and the 
high-energy description in terms of quark final states match onto each other very well. 

In fig. [6] (d), we find that, unsurprisingly, all limits improve as the dark sector gauge coupling 
a A' increases because the probability for A' radiation grows. Note, however, that in the very large 
a a' region, our perturbative treatment of the dark photon shower is expected to break down. 

Concerning the dependence of the sensitivity on the Z' mass (fig. [T] (e)), we find that at mz' < 
few TeV, the sensitivity of all searches increases with mz’- This is because a heavier Z' decays 
to more energetic x particles, which radiate more dark photons. Eventually, of course, on-shell 
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FIG. 8: 95% CL constraints on the dark photon mass rriA' and the kinetic mixing parameter e, with all other 
model parameters fixed at the first (second) set of benchmark values from table [I] in the left hand (right 
hand) panel. We show exclusion limits from the ATLAS search for prompt lepton jets in 5 fb -1 of 7 TeV 
data [15] (blue shaded region) and from their displaced lepton jet search in 20.3 fb^ 1 of 8 TeV data pi)] (red 
shaded region), as well as projected sensitivities for 100 fb -1 of 13 TeV data (blue/red unshaded regions). 
Black stars correspond to the nominal values of and e from table |T] The lighter colored region around 
m A < = 2 GeV corresponds to the transition region between the analysis in terms of hadron final states and 
the analysis in terms of quark final states and is based on interpolation. We also show the existing 90% 
CL exclusion limits from the electron and muon anomalous magnetic moment l66ll68| . HADES [69j, KLOE 
2013 [70] and 2014 [71], the test run results from APEX [75], BaBar 2009 [75] and 2014 [71], beam dump 
experiments E137, E141, and E774 [75H77] , A1 [75], Orsay [70], U70 [50], CHARM [5T], LSND [55], as well 
as constraints from astrophysical observations [55; [M] and 7r° decays [S3]. 


production of the Z l becomes impossible and the sensitivity decreases again. 

Finally, fig. [b] (f) shows that the sensitivity goes down when \ is made heavier because less 
boosted x particles radiate less. 

Let us now turn to a comparison of the two benchmark models from table [I] and compare figs. [6] 
and [7J The general features of the two figures are very similar, but in general benchmark point B 
is more easily detectable in prompt lepton jet searches due to the smaller cr. Benchmark point B 
offers somewhat better sensitivity than benchmark point A also because the smaller values of m x 
and mA' imply that the average number of A' radiated in each event is larger. This effect is 
especially pronounced in the displaced searches excluding type 2-2 events, which are limited by 
backgrounds. Note that in fig. [7] (c) the region with > 2 m x is not considered because it would 
lead to a very large branching ratio for A' —> xXi making the A' invisible to the detector. 

To put our results in the context of other constraints on dark photons, we show in fig. [8] the dark 
photon parameter space spanned by e and mA'- For both of our benchmark points, we compare 
the limits we derived from the prompt and displaced ATLAS lepton jet searches (blue shaded/red 
shaded) and the predictions for 13 TeV (blue/red unshaded) to the exclusion limits from various 
low energy experiments. For the 8 TeV displaced analysis, we decide for each parameter space 
point individually whether or not to exclude type 2-2 events (events with two calorimeter lepton 
jets), depending on which analysis leads to the better expected limit. 

We see that the parameter region probed by LHC lepton jet searches is complementary to the 
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region probed by low energy searches. It is important to keep in mind, though, that low energy 
experiments probe any dark dark photon model, while our analysis is sensitive only to scenarios 
that, in addition, feature a light DM particle that can be pair-produced at the LHC. We see that 
benchmark point A is already excluded by the 8 TeV ATLAS search for displaced lepton jets, while 
benchmark point B can be probed by both the prompt and the displaced search at 13 TeV. The 
exclusion regions from the displaced searches move to smaller values of e at larger rriA' because the 
A' width increases and the lab frame decay length decreases with increasing rriA' due to the smaller 
boost factors. Both effects need to be compensated by a decrease in e to avoid A’ decays in the 
inner detector. The impact of the various A' decay channels on the sensitivity is again reflected 
by spikes, dips and kinks in the excluded regions. For instance, note that the prompt search at 
7 TeV, which includes only muonic lepton jets, is insensitive close to the hadronically decaying p, 
uj and 0 resonances. At 13 TeV, the gap is closed thanks to the much better statistics. The sharp 
edge at low A' mass for the 7 TeV and 13 TeV regions corresponds to the opening of the decay 
channel A' — > fi + n~ at = ! 2m Remember that, in our 7 TeV and 13 TeV analysis, we could 
only include muonic lepton jets, therefore, our analysis is insensitive below the muon threshold. A 
full experimental search would of course extend also into this region, like the 8 TeV ATLAS search 
for displaced lepton jets has done. 


V. CONCLUSIONS 

In summary, we have studied the possibility of revealing the properties of a dark sector of particle 
physics by using final state radiation from dark matter produced at the LHC. The characteristic 
experimental signature of this process is a pair of lepton jets. While our conclusions apply to a 
very large class of extended dark sector models, we have worked in the framework of a toy model 
where fermionic dark matter particles \ are charged under a new dark sector gauge group U(l)' 
and coupled to the SM through a heavy mediator. When x particles are pair-produced at the 
LHC via this heavy mediator, they can radiate several light 17(1)' gauge bosons A 1 (dark photons) 
which subsequently decay to light SM particles (electrons, muons, mesons) through a small kinetic 
mixing e with the SM photon. Due to the required smallness of the kinetic mixing, the A' decay 
length can be macroscopic. We emphasize that our results for this toy model are easily generalized 
to any model with light dark sector particles charged under a new gauge interaction. We also 
remind the reader that, in order to account for all of the DM in the Universe, x must be produced 
non-thermally, as for instance in asymmetric DM scenarios. 

We have first developed two analytic treatments of dark photon radiation: in the first one, we 
use recursive expressions for the A' and x energy distributions, while in the second one, we compute 
the moments of these distributions fully analytically and then apply an inverse Mellin transform 
to obtain the distributions themselves. We have compared our analytic calculations with Monte 
Carlo simulations in Pythia, finding excellent agreement. 

In the second part of the paper, we have extended these Monte Carlo simulations by a simplified 
description of the ATLAS detector that allows us to recast the existing ATLAS searches for prompt 
and displaced lepton jets into powerful limits on the DM pair production cross section and the dark 
photon parameters in radiating DM models. Our limits on the XX production cross section range 
down to 0(10fb), depending on the A' mass and lifetime, the U(l)' gauge coupling, and the mass 
of X- Regarding the A' properties, we find that LHC searches for radiating DM can probe a region 
in niA '—e space that has been inaccessible to low energy dark photon searches so far, namely in 
the parameter range 10 - ' < e < 10~ 3 for m^/ in the MeV-GeV range. Of course, these limits are 
subject to the condition that DM particles with large coupling to dark photons (a a 1 ~ 0(0.1)) can 
be pair-produced at the LHC in significant numbers. 

Looking into the future, we have also shown that LHC limits will significantly improve in Run 2 
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at 13 TeV center of mass energy. It is worth emphasizing here that our simplified 13 TeV analyses 
are still a far cry from what a full experimental search can achieve. Most importantly, with a more 
detailed detector simulation and data-driven background estimation methods, it will be possible 
to include not only muonic lepton jets, but also A' decays to electrons and hadrons. 

To conclude, we have shown that lepton jets are an interesting and powerful tool to elucidate the 
dynamics of the dark matter sector. Together with searches for hadronic jets with non-standard 
properties, they will be essential in the hunt for dark matter at Run 2 of the LHC and might well 
be the first to find a positive signal. 
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